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Q-r ABSTRACT 

6 : 

;_h Aims. We combine near-to-mid-IR Spitzer data with shorter wavelength observations (optical to X-rays) to get insights on the properties of a 
^ sample of luminous, obscured Active Galactic Nuclei (AGN). We aim at modeling their broad-band Spectral Energy Distributions (SEDs) in 
C3 order to estimate the main parameters related to the dusty torus which is assumed responsible for the reprocessed IR emission. Our final goal 
1 1 ■ is to estimate the intrinsic nuclear luminosities and the Eddington ratios for our luminous obscured AGN. 

( Methods. The sample comprises 16 obscured high-redshift (0.9^z^2.1) X-ray luminous quasars (L2-io* f v~10 44 erg s~') selected from the 
, HELLAS2XMM survey in the 2-10 keV band. The optical-IR SEDs are described by a multi-component model including a stellar component 
' to account for the optical and near-IR emission, an AGN component which dominates in the mid-IR (mainly emission from a dusty torus heated 
\0 by nuclear radiation) and a starburst to reproduce the far-IR bump. A radiative transfer code to compute the spectrum and intensity of dust 
. reprocessed emission was extensively tested against our multiwavelength data. While the torus parameters and the BH accretion luminosities 
' are a direct output of the SED-fitting procedure, the BH masses are estimated indirectly, by means of the local M^/^-Mbu relation. 
, Results. The majority (~80%) of the sources show moderate optical depth (jgj IIm <S) and the derived column densities N# are consistent with 
■ the X-ray inferred values (10 22 1 N H 1 3xl0 23 cirr 2 ) for most of the objects, confirming that the sources are moderately obscured Compton-thin 
AGN. Accretion luminosities in the range 5x 10 44 1 L bol 1 4x 10 46 erg s are inferred from the multiwavelength fitting procedure. We compare 
t-H ■ model luminosities with those obtained by integrating the observed SED, finding that the latter are lower by a factor of ~2 in the median. The 
JL* ■ discrepancy can be as high as an order of magnitude for models with high optical depth (t 91p „, = 10). The ratio between the luminosities 
obtained by the fitting procedure and from the observed SED suggest that, at least for Type 2 AGN, observed bolometric luminosities are 
^ . likely to underestimate intrinsic ones and the effect is more severe for highly obscured sources. Bolometric corrections from the hard X-ray 
5^ ' band are computed and have a median value of k 2 _iofa.y~20. The obscured AGN in our sample are characterized by relatively low Eddington 
ratios (median /t £( M~0.08) . On average, they are consistent with the Eddington ratio increasing at increasing bolometric correction (e.g. 
IVasudevan & Fabian 200jJ. 
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1. Introduction 

A robust determination of Active Galactic Nuclei (AGN) 
Spectral Energy Distributions (SEDs) is of paramount impor- 
tance to better understand the accretion processes onto su- 
permassive black holes (SMBHs) and their cosmological evo- 
lution. According to our present knowledge, the bulk of ac- 
cretion luminosity is emitted in the optical-UV range with a 
quasi-thermal spectrum originating in an optically thick, geo- 



metrically thin, accretion disk. Electrons with temperatures of 
the order of a few hundreds of keV form a hot corona which 
upscatters disk photons to X-ray energies with a power law 
spectrum and an exponential high- energy cut-off correspond - 
ing to the electron temperature (e.g. jHaardt & Maraschill991 ). 
Dusty material, possibly with a toroidal geometry, intercepts a 
fraction of the primary continuum which depends on the cover- 
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ing factor. The absorbed energy is re-emitted in the near to far 
infrared with a grey-body spectral shape. 

The SED of optically bright, unobscur ed QSO is relativel y 
well known. After the seminal work of Elvis et al.1 dl9 94). 
fairly a ccurate measurements were published bv lRichards et al.l 
(120061) using Sloan Digital Sky Survey (SDSS) data. The emis- 
sion is characterized by a double bump. In a vF v diagram, the 
optical-UV spectrum rises steeply towards the shortest acces- 
sible wavelengths. It is commonly referred to as the Big Blue 
Bump and thought to be the accretion disk signature. The IR 
bump is weaker and likely due to a dusty torus seen almost 
face-on. The ratio between UV (at 2500 A) and X-ray lumi- 
nosity (at 2 keV), parameterized by the slope or J3 of the power 
law connecting the rest-fram e luminosities, increases with in- 
creasi ng UV luminosity (e.g., Vignali et al. 20031; Isteffen et al 
20061) . 



The average S EDs for rad i o-loud and radio-quiet Type 1 
AGN presented in Elvis et al.l (119941) allow the estimation of 
bolometric corrections, which are a key parameter to deter- 
mining the bolometric luminosity from observations at a given 
frequency, and the Eddington ratio, once the SMBH mass 
is kno wn. By includ i ng the a ox vs. UV luminosity depen- 
dence, iMarconi et a D d2004 and iHopkins et all (|2006) com- 
puted luminosity-dependent bolometric corrections and adopt 
them to estimate the local SMBH mass density from the ob- 
served X-ray luminosity functions. 

The luminosity depen dence of bolometric corre ctions 
was recently questioned by IVasudevan & Fabian] d2009l) . who 
pointed out the importance of simultaneous optical-UV and 
X-ray observations and reddening corrections in the UV. They 
suggest that the bolometric correction correlates with the 
Eddington ratio rather than with the bolometric luminosity. 
Their observational resul ts align with th e predictions of ac- 
cretion disk models (e.g. JWitt et alj[l997l) where higher blue- 
bump to X-ray ratios for sources with higher Eddington ratios 
are expected. 

While important progress has been made towards a bet- 
ter determination of Type 1 AGN SEDs, our knowledge of 
Type 2 broad-band spectra is much more limited despite the 
fact that most of the accretion-driven energy density in the 
Universe is expected t o occur in obscured AGN (e.g., [Fabian 
1999llGilli etaDl2007 and references therein). Therefore, a ro- 
bust estimate of their bolometric luminosity is extremely im- 
portant to properly address the issue of SMBH evolution over 
cosmic time. 

Nuclear accretion luminosity in Type 2 AGN is very faint 
in the optical-UV and soft X-rays. Moreover, the host galaxy 
stellar light often dominates in the optical making it difficult to 
disentangle nuclear emission from starlight. Infrared emission 
is only marginally affected by dust obscuration and has proved 
to be a powerful indicator of dust obscured AGN. In particular, 
the thermally reprocess ed nuclear emis sion of obscured Type 2 
AGN is expected (e.g jFritz et al . 2006, hereafter F06, and ref- 
erences within) to peak around a few tens of micron. 

Mid-IR (MIR) observations, and especially those obtained 
in the last few years with the Spitzer satellite, are extremely 
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efficient in the study of obs c ured AGN (e.g. | Rig bv et al 
2005 : iMartfnez-Sansigre et al.l 120051: IWeedman et al] 120061: 
Fior e et al. l2008h . In a previous paper ( Pozzi et ail 2007h we 



presented the first analysis of the mid-IR data of a Spitzer 
program devoted to a systematic study of the broad-band 
properties of X-ra y selected, luminous obscured quasars. In 
Pozzi et alJ (120071) . the SEDs were re produced by means of 
SED templates from lSilva etail d2004 . 

Here we present the observational data for our final sam- 
ple of 16 obscured quasars and the detailed modeling of 
their broad-band SED using a more complete multi-component 
model, with goodness of fit estimated via a% 2 analysis. 

The outline of the paper is as follows: in §2 the X-ray se- 
lected quasar sample is presented, along with all the available 
multi-band (optical, near-IR (NIR) and sub-mm) and spectro- 
scopic follow-up. The Spitzer data are presented in §3, along 
with data reduction and analysis. In §4, the complete multi- 
component model is described. In §5, the best-fitting solutions 
are discussed, while in §6 we focus on the black hole physical 
properties that can be constrained from the best-fitting proce- 
dure. Finally, the main results are summarized in § 7. 

Hereafter, we adopt the concordance cosmology (Hp = 
70 km sec 1 Mpc 1 , Q m =0.3 and Q A =0.7. lspergel et all2003l) . 
Magnitudes are expressed in the Vega system. 



2. The sample 

The sample presented in this work comprises 16 X-ray 
obscured quasars detected in the HELLAS2XMM survey 
(Baldi et al.l 120021) and observed by Spitzer in 2006. The 
HELLAS2XMM survey is a shallow, large-area hard X-ray sur- 
vey (S2-\okeV > 10~ 14 ergcirT 2 s _1 ) over a final area of 
1.4 deg 2 . The catalogue comprises 232 X-ray sources; ~92% of 
the sample is optically identified down to R~25, while ~70 % 



of the sources have a spectroscopic classificati on dFiore et al 
2003, Maiolino 'eTaIll2006llCocchia et al.ll2007h . 



The 16 sources were selected from the original survey in 
order to include the most luminous obscured quasars. The se- 
lection was primarily based on the X-ray-to-optical flux ratio 
(hereafter X/oJH, which has been proved to be an e fficient way 



to sel ect high-redshift (z; 1), obscured quasars (see lFiore et al 
2003). All but one of the sources were selected to have X/O 



greater than 1 (see Fig. Q]and Table [TJ; the only exception is 
GD 158#19 (X/O-0.63) w hich was included i n the sample for 
its peculiar properties (see Vignali et al. 20091 hereafter V09). 
We note, however, that not all of the HELLAS2XMM sources 
matching this selection criterion are present in this work. 

The selected sources are relatively faint in the optical band, 
with an R-band magnitude in the range 21.8-25.1 (the brightest 
object being the peculiar source GD 158#19). By combining 
optical photometry with deep K s -band photometry (obtained 
with the Infrared Spectrometer And Array Camera, ISAAC, 
mounted on the ESO-VLT1 Telescope), almost all the objects 
are found to be Extremely Red Sources (EROs, R-K s >5). The 
link between high X/O ratios and optical-to-near colours was 



X/O 



is defined as log . We used fitiO) 
1.74X10-' erg cnrV 1 A" 1 and AA R = 2200 A. IZomberid dl990h . 
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Fig. 1. R-band magnitudes vs. har d X-ray (2-10 keV) fl ux for 
the full HELLAS2XMM sample dCocchia et al.ll2007b . Blue 
triangles represent the sources included in the present anal- 
ysis: blue triangles inside red symbols represent the sources 
with spectroscopic redshifts. Other symbols: empty red squares 
= sources spectroscopically classified as Type 2 AGN; empty 
circles = sources spectroscopically classified as non-Type 2 
AGN (Type 1 AGN, emission-line galaxies, early-type galaxies 
and groups/clusters of galaxies); crosses: objects not observed 
spectroscopically; upward arrows = lower limits. The dashed 
lines represent the loci of constant X/O ratio (X/0=+l). 



studied by e.g., Bru sa et al.1 (12005). Considering different X-ray 
surveys at different depths, they find a clear trend: the higher 
X/O, the redder the source. 



In iMignoli et al.1 ([2004), detailed K,-band morphological 
studies were presented for 8 objects of the sample selected from 
those sources with the more extreme R-K, colours: the majority 
of the sources (6 over 8) have an extended K s -band morphol- 
ogy, consistent with an elliptical-type profile, without evidence 
for a nuclear point-like source (which would be expected to 
trace the X-ray AGN). This suggests that the nuclear emission 
is diluted and hidden by the host galaxy up to at least 2.2pm. 
In Table [TJ the R and K s -band magnitudes are reported. 

For 11 sources, spectroscopic in formation i s avai lable 
thanks to optical (8 sources, see Fiore et al. [2003 and 



Cocchia et al 



Maiolino et a 



2007) and near-IR spectroscopy (3 sources, see 



2006 and Sarria et al., in preparation). All but 



one spectra are typical of optically obscured AGN, thus con- 
firming the X-ray classification (i.e., Type 2). The only excep- 
tion is source PKS 0537#91, with emission line ratios typi- 
cal of an HII region (Sarria et al. in preparation). One source, 
Abell 2690#29, shows the typical rest-frame spectrum of a 
high-redshift dust-redden ed quasar, with a bro ad Ha line and 
a Type 1.9 classification dMaiolino et al.ll2006l) . Spectroscopic 
redshift z are in the range 0.9-2.08 and are reported in TableQ] 

For the sources witho ut redshift (5 out of 16), a photometric 
redshift was estimated in lPozzi et al. ( 20071) . where the Spitzer 
data reduction and a preliminary SED analysis were presented. 



Finally, sub-mm observations were performed for four 
sources in October-November 2004. Only one object, 
GD 158#19 (z=1.957), was detected, while for the others up- 
per limits were gained. Because of its broad-band coverage 
(up to 850/mi) with good-quality photometric data, the source 
GD 158#19 was studied in a dedicated work (see V09). 

In Table [TJ we report source name, 2-10 keV flux, R, and 
K s photometry, X/O ratio, the redshift z, the column densities 
N// and the absorption-corrected (2-10 keV) X-ray luminosi- 
ties of the sample. The source table order reflects the Spitzer 
observation strategy (see Sec 3.). The sub-mm flux densities 
are reported in Table[2]along with the fluxes obtained in the IR 
bands (from 3.6pm up to 160/im) with the Spitzer satellite (see 
§0. 

Almost all the sources have column densities N# in the 
range 10 22 -10 23 ' 4 crrT 2 and 2-10 keV rest-frame luminosi- 
ties in the range 10 43 8 -10 44 7 erg s _1 , placing them among the 
Type 2 quasar population. 



3. The Spitzer data 

The targets were observed by Spitzer in 2006 with both IRAC 
and MIPS instruments in photometry mode. All the sources 
were observed with the same total integration time in the IRAC 
bands (480s), while different strategies were followed in the 
MIPS bands taking into account the different optical-NIR prop- 
erties. While all the objects were observed at 24pm for a total 
integration time of 1400s, only a sub-sample of sources with 
relatively bright R-band magnitudes (R<24) and spectroscopic 
redshifts were observed at longer wavelengths, with integra- 
tion times of 300s and 600s at 70 and 160/im, respectively. 
PKS 0537#91 (see Table [TJ was not observed at 70 and 160 
pm since it had no redshift at the epoch of the Spitzer observa- 
tions. 

The reduction method is described in detail in Pozz i et al 



( 2007) and in V09 and it is briefly summarized here. The IRAC 
flux densities of the sources were measured from the post- 
basic calibrated data (post-BCD) images in the Spitzer archive. 
Aperture fluxes were measured on the background subtracted 
maps within a 2.45 aperture radius using aperture corrections 
of 1.21, 1.23,1.38 and 1.58 for the four IRAC bands (following 
the IRAC Data Handbook). For the MIPS bands, we started 
the analysis from the basic calibrated data (BCD) at 24pm 
and from the median high-pass filtered BCD (fBCD) at 70 and 
160/im, as suggested for faint sources. At 24pm, the BCD were 
corrected for a r esidual flat fielding dependent on the scan mir- 
ror position (see Faddaetal . 2006; Pozzi et al ] l2007l) . We con- 
structed then our own mosaic s using the SSC MOPEX software 
(Ma kovoz & Marle au 2005). Aperture fluxes were measured 
within a 7 aperture radius for the 24pm band and 16 aper- 
ture radius for the 70 and 160/im band. The aperture corrections 
used were 1.61, 2.07 and 4.1, respectively (see the MIPS Data 
Handbook). A small aperture radius was used at longer wave- 
lengths (at 160//m the adopted radius is comparable to half of 
the PSF FWHM) to exclude the contamination of nearby far- 
infrared sources (see V09). At 24pm, thanks to a better PSF 
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Table 1. Properties of our luminous obscured quasars 



Complete Name 


Name 


2-10 keV flux 


R 


K s 


X/O 


z 


N H 


^2-10 keV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 

1 n c9.4 


(9) 



HELLAS2XMM 053920.4-283721 PKS0537#lla 4.19 23.40 18.25 1.48 

HELLAS2XMM 053917.1-283819 PKS 0537*164 1.50 23.60 19.02 1.12 

HELLAS2XMM 053851.3-283949 PKS 0537*123 2.97 23.10 17.94 1.21 

HELLAS2XMM 003413.8-115559 GD 158*62 3.59 23.30 21.83 1.38 

HELLAS2XMM 003357.2-120039 GD 158*19 2.43 21.80 - 0.61 

HELLAS2XMM 204428.7-105629 Mrk 509*01 2.10 23.85 17.88 1.36 

HELLAS2XMM 204349.7-103243 Mrk 509*13 3.10 23.99 18.79 1.59 

HELLAS2XMM 235956.6-251019 Abell 2690*75 3.30 24.60 18.33 1.85 

HELLAS2XMM 031343.5-765426 PKS 0312*36 1.90 24.70 19.13 1.66 

HELLAS2XMM 054021.1-285037 PKS 0537*91 4.2 23.70 18.99 1.60 

HELLAS2XMM 053945.2-284910 PKS 0537*54 2.1 25.10 18.91 1.86 

HELLAS2XMM 053911.4-283717 PKS 0537*111 2.1 24.50 17.64 1.62 

HELLAS2XMM 0001 11.6-251202 Abell 2690*29 2.8 25.10 17.67 1.99 

HELLAS2XMM 031018.9-765957 PKS 0312*45 2.8 24.40 18.62 1.70 

HELLAS2XMM 005030.7-520046 BPM 16274*69 2.27 24.08 17.87 1.48 



0.981 
1.824 
1.153 
1.568 
1.957 
1.049 
1.261 

i 3+0.30 t 
1- -0.20 
n q+0.05 t 
u ' 7 -0.15 

1.538 
>1.3 f 

i 7+0.20 t 
L ' -0.10 

2.08 
1 85 +a2 ° f 
1.35 



l "- / 4.8 

1 3 15 



6.6-" 
26.3£I 
7 3 1L7 

<1.1 

2.5« 
1 L2 

A.S Q'020 



9 l' 2 - 4 
2-5|;S 



1.9 
3.1 
2.0 
4.9 
6.3 
1.5 
2.6 
3.2 
0.7 
8.1 
2.0 
1.7 
8.4 
6.2 
2.4 



(1) Source complete name; 

(2) source abbreviated name (adopted throughout the paper); 

(3) 2-10 keV X-ray fluxes in units of 10~ 14 erg cm -2 s" 1 from lPerola et al.1 [2004 ) (with the exception of BPM 16274*69, from Lanzuisi et al., 
in preparation); 

(4) R-band magnitude from lFiore et al.l j2003n (with the exception of source BPM 16274*6 9. fromlCocchia et alj2007h ; 

(5) K s -band magnitudes. For a sub-sample of sources, the Kj-band analysis can be found in lMignolietalJ ( l2004l) ; 

(6) X-ray-to-optical flux ratio; 

(7) source redshifts. Spectroscopic redshift from optical spectros copy dFiore et al.ll2003h and, for three sources, from near-IR spectroscopy 
(Abell 2690*29 and BPM 16274*69 from iMaiolino et afll2006l : PKS 0537*91 from Sarria et al., in preparation); photometric redshifts 
( T symbol), and corresponding lcr errors, from pTzzTeTalT i 20071) ; 



(8) column densities in source rest-frame in units of 10 cm measured from X-ray spectral fitting (see Perola et al . 2004 and Lanzuisi et al., 
in preparation). For the sources with photometric redshifts, values are taken from lPozzi e t al. (2007). For PKS 0537*164 the X -ray spectral 
fitting is prevented by the few X-ray counts; for Mrk 509*01, only an upper limit was derived dPerola et al.ll2004l) ; for PKS 0537*54, the 
column density is 9.3xl0 22 cm" 2 for z = 1.3. Galactic absorption column densities adopted in the spectral fitting are: 8xl0 20 cirT 2 for the 
field PKS 0312-77; 2xl0 20 cnT 



2.5xl0 20 cm' 2 for the field GD 158-100 (see lStark etai]|l992h ; 

(9) rest-frame 2-10 keV absorption-corrected X-ray luminosity in units of 10 44 erg s _1 from lPerola et alj |2004) and from lPozzi et al.l d2007h 
for sources with photometric redshifts. Luminosities are computed using H =70 Km s _1 Mpc~', Cl m = 0.3 and Q A = 0.7. 



sampling, two sources (PKS 03 12#36 and Abell 2690#29) were 
deblended using a PSF deconvolution analysis. 

All 16 sources were clearly detected in the IRAC bands. 
At 24/im, 14 sources (out of 16) were detected above the 5cr 
level and span almost two orders of magnitude in flux, from 
~7000 jj]y down to the faintest source, close to the 5<x detec- 
tion level (~100 /uJy). For the two sources without detection, 
an upper limit (3<x) was estimated from the average noise of 
the map, derived by making multiple aperture measurements at 
random locations throughout the residual mosaic after source 
extraction. The typical average noise (lcr) is 20 fiiy. At 70 and 
160^m, as said before, only the brightest R-band sources were 
observed; among them, only the two most luminous (in the op- 
tical band) were detected, PKS 0537#43 and GD 158#19, the 
latter being the source described in V09. For the remaining 6 
sources, an upp er limit (3<x) was estimated from the residual 
mosaic (see also iFraver et al.ll2006l) after source extraction, as 
done at 24yum. The typical average noise (lcr) is 1.2 mJy and 



8 mJy at 70 and 160//m, respectively (c onsistent with the results 
obtained in the COSMOS field from lFraver et al.ll2009l taking 
into account the different integration time). 

Table[2]reports the target flux densities provided by Spitzer. 
To compute uncertainties, the noise map was added in quadra- 
ture to the systematic uncertainties, assumed to be 10 per cent 
in the IRAC and MIPS 24yum bands and 15 per cent at 70 and 
160/Lan (see IRAC and MIPS Data Handbook). 



4. Modeling the spectral energy distribution 

The observed optical-to-MIR (or FIR/sub-mm) SEDs can be 
modelled as the sum of three distinct components: a stellar 
component, which emits most of its power in the optical/NIR, 
an AGN component, whose emission peaks in the MIR for ob- 
scured quasars, and is due to hot dust heated by UV/optical ra- 
diation from gas accreting onto the central SMBH, and a star- 
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Source name 


3.6pm 


4.5pm 

L"Jy] 


5.8pm 

[j"Jy] 


8.0pm 


24pm 

0"Jy] 


70pm 
[mJy] 


160pm 
[mJy] 


450pm 
[mJy] 


850pm 
[mJy] 


PKS 0537#43 


423 ± 42 


918 i 


:92 


1750+175 


3202 + 320 


6879 ±688 


27.0±4.5 


39.6± 7.8 






PKS 0537#lla 


57 


± 6 


69 


± 7 


95± 


10 


168 ± 


17 


517 +57 


<3.6 


<24 






PKS 0537*164 


21 


± 2 


18 


± 2 


19+4 


23 : 


t 5 


<60 


<3.6 


<24 






PKS 0537*123 


84 


±9 


96 d 


: 10 


113+ 


13 


165 + 


17 


747 ±81 


<3.6 


<24 






LjIJ iJotfuZ 


60 


±6 


91 


+ 9 


155± 


17 


339 ± 


35 


lOUo ±104 


< J.O 


<Z4 


< /4 


^6. 1 
<0. 1 


GD 158*19 


226 i 


23 


387 h 


:39 


756± 


76 


1547 ± 155 


5326 ±534 


27.0±4.5 


37.7± 7.8 


<94.3 


8.6±2.1 


Mrk 509*01 


63 


± 7 


51 


+ 6 


60± 


11 


59 + 


10 


<60 


<3.6 


<24 


<82 


<6.6 


Mrk 509*13 


51 


± 6 


69 


± 7 


110+ 


11 


215 + 


17 


866 ± 94 


<3.6 


<24 


<75 


<6.2 


Abell 2690*75 


51 


± 5 


56 


± 6 


89± 


11 


139 + 


15 


565 ± 62 










PKS 0312*36 


41 


± 4 


44 


+ 5 


40+ 8 


710 : 


t9 


236 ± 30 










PKS 0537*91 


28 


± 4 


35 


±4 


42 


t 8 


800 ± 


10 


301 ± 40 










PKS 0537*54 


31 


± 4 


35 


± 4 


50+ 


10 


470 : 


t 8 


279 ± 45 










PKS 0537*111 


88 


±9 


75 


± 8 


41 


t 6 


460: 


t 7 


148 ± 28 










Abell 2690*29 


141 ± 14 


185 h 


: 19 


260± 


27 


371 + 


38 


1012 ± 106 










PKS 0312*45 


50 


±6 


62 


± 7 


69+ 


10 


780 ± 


10 


249 ± 35 










BPM 16274*69 


86 


±9 


92 


+ 9 


97± 


11 


120 ± 


13 


286 ± 34 











The upper limits are given at the 3cr confidence level. 
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Fig. 2. x 2 distribution. The red hatched and the empty distribu- 
tions represent the best-fitting solutions and all the solutions at 
lcr, respectively. PKS 0537#43 is not reported given the high 
X 1 value (Table [3}. 



burst component, which represents the major contribution to 
the FIR spectrum. 

In this work, we considered all the three components (see 
Sec. 4.1, 4.2 and Hatziminaoalou et al. 2008). Since the focus 



of the paper is on the AGN contribution to the SED, we discuss 
in some more detail the hot dust modeling and its uncertainties. 
The hot dust emission in AGN is reproduced using the F06 



ent authors (e.g., Pier & Krolik 


19921 Granato & Danese!l994: 


Efstathiou & Rowan-Robinson 


1995), where the IR emission 



in AGN originates from dusty gas around the SMBH with a 
smooth distribution. The dust grains are heated by high-energy 



photons coming from the accretion disk, and their thermal and 
scattering re-emission, mostly at IR wavelengths, is computed 
by means of the radiative transfer equations. For what concerns 
the dust distribution geometry, different possibilities (i.e. "clas- 
sical" torus shape, tapered or flared disk), are explored in the 
literature. 

More recently, models considering a clumpy distribution 
for the dust have been developed (e.g., iNenkova et al] l2002t 



Honig et al. 2006; N enkova et alj 12008). These models sue 



cessfully explain many recent observations in the mid-IR, such 
as the strength of absorp tion and emission f eatures at 9.7 pm 
and the X-ray variability (Risa liti et al ] l2002h . 

A further possibility to the t orus models described a bove 
are the disk-wind models (see Elitzur & Shlosman 2006 and 
references therein), involving a completely different approach. 
The dusty clouds, responsible for the obscuration, are part of an 
hydro-magnetic wind coming from the accretion disk. On the 
one hand, these models are potentially capable of explaining 
disparate phenomena in AGN (from broad emission to absorp- 
tion lines and obscuration), providing an hydro-dynamicaljus- 
tification for the persistence of the clouds around the SMBH. 
On the other hand, a parameterization which takes into account 
the observational constraints on the clumpy obscuration, hence 
supplying a g rid of synthetic IR SEDs, does not exist yet (see 
lElitzurll2008l) . 

A detailed comparison betwee n smooth and clumpy dust 
distrib ution models is discussed by Dullemond & van Bemmell 
(2005J). They conclud e th at both models yield similar SEDs 
(see also lElitzurll2008l and Nenkova et al . 2008). The main dif- 
ference is in the strength of the silicate feature observed in ab- 
sorption in objects seen edge-on, which is, on average, weaker 
for clumpy models with the same global torus parameters. In 
clumpy models, in fact, clouds at different distances from the 
central source can be intercepted by the line of sight, including 
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the innermost clouds, where the silicate feature is in emission 
given the higher temperature of the dust grains. 

A systematic comparison of the two model predictions is 
beyond the scope of the present paper and should be performed 
on high-quality IR data (i.e. a Spitzer IRS spectroscopic sam- 
pie). 

Notwithstanding these limitations, with the present work 
we aim at extracting the maximum information using the avail- 
able photometric data. The F06 model adopted in this work 
is one of the mo dels best tested a gainst both broad-band 
photometry (F06: iBerta et al.1 120041 iRodighiero et all l2007t 
Hatziminaoglou et al.ll2008l 120091: Kgol et al.ll2009t V09) and 
Spitzer mid-infrared spectra (F06). Moreover, the F06 model 
was the first able to reproduce the quasar mid-IR spectra, con- 
sidered a very important constraint to characterize the dust 
properties in AGN and probe the Unified Model. 

4.1. The AGN-torus component 

The F06 code assumes a smooth dust distribution around the 
central source consisting of a Galactic mixture of silicate and 
graphite grains. The presence of silicate dust grains is clear 
from the absorption feature at 9.7 fim seen in most of Type 2 
AGN. The graphite grains are, instead, responsible for the rapid 
decline of the emission at wavelengths shortwards of a few mi- 
crons, corresponding to a blackbody emission of about 1500 K, 
the sublimation temperature of these grains (see F06). 

The assumed dust geometry is a flared disk (see 
Efst athiou & Rowan-Robinsonll 19951) . that is a sphere with the 
polar cones removed. The internal radius of the dust distribu- 
tion is defined by the sublimation temperature of the dust itself. 
In order to simulate a more realistic shape for the dust distribu- 
tion, F06 assumes that the dust density can vary both with the 
radial and the angular coordinates: 



p(r,0) =po-r° e 



-y|cos(9)| 



(1) 



where 6 is the angle with respect to the equatorial plane. 

The dusty torus is heated by the emission of the inner ac- 
cretion disk, which represents the input energy to the radiative 
transfer code. The assumed spectrum for the accretion disk is 
defined in the 10~ 3 to 20/mi regime (from soft X-rays, i.e. 1.25 
keV, to mid-IR) and is parameterized by broken power laws in 
F06. The specific indices of the power laws are adapted from 
the Granato & Danese (1994) and the Nenkova (2002) mod- 
els and are consistent with the broad-ban d SEDs of a sample 
of Ty pe 1 AGN from the SDSS survey dHatziminaoglou et al. 
2008; see their Figs. 8 and 9). 

Along with the thermally re-processed light, the F06 pro- 
vides, as a function of the line-of-sight and optical depth, the 
fraction of the inner accretion disk light not intercepted by the 
torus and the scattered light. In the following, with AGN com- 
ponent we will refer to the sum of all the three contributions. 

4.2. The stellar and starburst components 

The stellar component is modelled as the sum of Simple Stellar 
Populations (SSP) models of different age, all assumed to have 



common (solar) metallicity. A Salpetei ( 19551) initial mass 
function (IMF) with mass in the range (0.15-120 M ) is as- 
sumed. The SSP spectra h ave been weighte d by a Schmidt-like 
law of star formation (see lBerta et al.ll2004l) : 



SF/?(0 = Zkzi X exp 1- Tg ~ 1 



T g 



(2) 



where Tq is the age of the galaxy (i.e. of the oldest SSP) and 
r s f is the d uration of the burst in units o f the oldest SSP (see 
V09). As in lHatziminaoglou et al.l (12008b and V09, a common 
value of extinction is applied to stars at all ages, adopting t he 



Cardelli et al. 1989). 



Hop 

aLL 



extinction law of our Galaxy (Ry = 3.1; 

In order to keep the number of free parameters as low as 
possible, emission from cold dust, which dominates the bolo- 
metric emission at wavelengths longer than 30/mi rest-frame, 
is included only when far-IR/sub-mm data allow us to con- 
strain that part of the SED (two sources of the sample). For 
the same reason, additional comp onents, such as the c old ab- 
sorber detached from the torus (i.e. Pollett a et al . 2008), which 
might improve the fit but would increase the complexity of the 
overall modeling, is not included. To reproduce the starburst 
component, a set of semi-empirical models of well known and 
studied starburst galaxies is used, similarly to V09. 

4.3. SED fitting procedure 

The quality of the fitting solutions is measured using 
a standard y 1 minimizat ion technique (as also done in 
Hatziminaoglou et al ■Il2008l) . where the observed values are the 
photometric flux densities (from optical-to-MIR/FIR) and the 
model values are the "synthetic" flux densities obtained by 
convolving the sum of stars, AGN and starburst components 
through the filter response curves. 

Before starting the general fitting procedure, we test which 
torus parameters mainly influence the global model SED and 
are more sensitive to our data sets. Parameters which are not 
constrained by our data were then frozen. 

The F06 torus model is described by six parameters: the 
ratio R m axlRmin between the outer and the inner radius of the 
torus (the inner radius being defined by the sublimation tem- 
perature of the dust grains); the torus full opening angle 0; the 
optical depth r at 9.7/vm (T9.7); the line of sight 9 with respect 
to the equatorial plane, and two parameters, y and a, describing 
the law for the spatial distribution of the dust and gas density p 
inside the torus (Eq.Q}. 

In our approach, we let free to vary, inside the pre- 
constructed grid of torus models the following parameters: the 
torus full opening angle ©, the optical depth 79,7 and the param- 
eter a describing the radial dependence of the density. We fix 
RmaxlRmin-^, which translates into compact tori of a few tens 
of parsecs (given that /?„„•„ is directly connected to the sublima- 
tion temperature and to the accretion luminosity of the central 
BH). Recent high-resolution IR observations support a com- 
pact dust distribution in nearby luminous AGN. Using the in- 
terferometry at VLTI in the 8-1 3 ^urn band, a tor us of size ~2-3 
pc was detected in NGC 1068 djaffe et al.ll2004 . Similar com- 
pact tori were also found in other local AG N, as Circinus and 
NGC 4151 (see the review bv lElitzurll2008l) . 
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Fig. 3. a. Observed-frame SEDs for 14 sources with data from the R-band to the 24/im (black dots) compared with the best-fit 
model obtained as the sum (solid black line) of a stellar component (red dotted line) and an AGN component (blue dashed line). 
The nuclear Kj-ban d upper limits (downward-pointing arrows) were derived from the morphological analysis carried out by 
Mignohetal](l2004 . 



Regarding the density distribution, we allow power-law 
profiles decreasing with the radius with different coefficients 
a. We do not allow any dependence on the distance from the 
equatorial plane by fixing y=0. As a result, different angles for 



the lines of sight 8 (with respect to the equatorial plane) give 
the same SED, once the torus is intercepted. 

Given the F06 grid of models, the discrete values, allowed 
for our free parameters, are: = [60, 100, 140°], 79,7 = 
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[0.1,0.3,0.6, 1,2, 3, 6, 10] and a = [-1, -0.5, 0], implying 72 
different torus SEDs. 

Considering the stellar component, we fix at z — 4 the red- 
shift for the formation of the oldest SSP, i.e., given the observed 
redshift of the sources, we consider galaxies with ages typical 
of early-type galaxies (I 1-2 Gyr). This assumption is justified 
by the observed R-K, colours and the brightness profiles typi- 
cal of early-type galaxies as obtained by a detaile d morpholog- 
ical analysis in the K s -band (Mi gnoli et al.1 2004). Concerning 
the star-formation history, we allow the T s f parameter of the 
Schmidt-like law and the value of the extinction E(B - V) to 
vary; the latter is a key parameter, along with the optical depth 
Tg .7 of the torus, in shaping the optical-NIR continuum. 

Overall, the SED-fitting procedure ends with 5 free pa- 
rameters (6 when also a starburst component is included). 
Since the problem is affected by some degeneracy, we consider, 
along with the best-fitting solutions, all the acceptable solutions 
within lcr confidence level (by considering, for each source, all 
the solutions with^ 2 -xti n = Ay 2 ! 5.89 or 7.04, when the cold 



starburst component is added, see lLampton et al ]fl97rl 



5. Results from SED fitting 

In Fig. [2] we show the^- 2 distribution. The hatched histogram 
represents the distribution of the best-fitting solutions for our 
16 targets, while the empty histogram shows the x 2 distri- 
bution of all the solutions within lcr, satisfying the criteria 
X 2 ~Xmi„ + Ax- Given our adopted grid for the fitted parameters, 
the total number of solutions at lcr (constructed by adding all 
the solutions at lcr of each object) is 137, 8 solutions on aver- 
age per source (including the best-fitting one). 

The two distributions do not show a significant difference, 
and this reinforces our choice of considering, in the analysis of 
the parameter space and degeneracy, all the solutions at lcr as 
a unique statistical sample. 

In terms of the absolute values of the x 2 , only 7 sources 
(out of 16) give a formally acceptable fit (P(x 2 > X 2 b s ) > 90%, 
see Table |3); the remaining sources have a best-fit model with 
large X 2 b s - While we use the x 2 to assign a relative goodness 
of different parameter combinations inside the parameter grid, 



we will not take the absolute probabilities at face value. Over- 
estimated x 2 are, in fact, a common problem of most SED fit- 
ting techniques, resulting from a combination of two different 
reasons: the limited grid of models (72 torus models with the 
adopted choice of parameters, see 94.3b with no uncertainties 
associated, and the photometri c measurements with o ften un- 
derestimated uncertainties (see iGruppioni et al. 2008 for a de- 
tailed description of this issue). 

In Fig.|3^,b the observed SEDs, from the R-band to the IR 
(or sub-mm), are reported with the best-fitting models over- 
plotted (solid line). All the sources need a host galaxy compo- 
nent (red dotted line) and an AGN one (blue long-dashed line). 
The stellar component dominates in the R and K, bands, while 
the nuclear one at 24/rni. In the IRAC bands, both components 
contribute, the fraction depending on the properties of the in- 
dividual sources. For PKS 0537.43 and G158.19, where data 
points at longer wavelengths are available, an additional star- 
burst component is needed (green dot-dashed line in Fig.|3p). 

In Fig. |4] the relative contributions of the thermal, di- 
rect and scattered light to the total AGN light are shown for 
two sources characterized by a low (t9.7=0.1, PKS 0537#123) 
and a high (t 9 . 7 =10, PKS 0537#111) optical depth, respec- 
tively. While for r 9 7=10, the AGN emission is dominated by 
the reprocessed emission in all the UV/optical/IR bands, for 
T9.7 =0. 1 the direct and scattered components account for the 
optical/UV AGN emission. Nevertheless, the contribution of 
the components mentioned above never exceeds the 20% of 
the observed flux in the R-band. For a sample of highly polar- 
ized red active galactic nuclei selected from the 2MASS sur- 
vey, a larger contribution of the scattered nuclear component to 
the optical and near-IR e mission was found (Cutri et al. 2002, 
Kuraszkiewicz et al]|2009h . The different result obtained from 
our analysis is probably due to the different degree of obscura- 
tion of the two samples: the 2MASS sample is characterized by 
column densities around 10 22 cirT 2 , while the present sample 
has a median column density of 7xl0 22 crrT 2 . 

The negligible contribution of the AGN component, rela- 
tive to the stellar one, at short wavelengths, is c onsistent with 
the upper limits to the AGN emission derived bv lMignoli et al 
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Fig. 4. Observed-frame spectral energy distribution as in Figs. 
3a,b for 2 sources characterized by extreme values of the 
optical depth: r 9 . 7/im =0.1 and 10 for PKS 0537.123 and 
PKS 0537_1 1 1 respectively. In this Figure, the 3 different com- 
ponents that contribute to the AGN emission are reported: di- 
rect nuclear light (long-dashed line), scattered light (dashed-dot 
line) and thermally re-emitted light (dashed-dot-dot-dot line). 



(2004) from the analysis of the K s -band images (shown as 
downward-pointing arrows in Figs.[3]and|4]l. 

In Table [3] the x 2 values (and the corresponding degree 
of freedom) of the best-fitting solutions are reported for each 



5.1. Torus parameters 

In the following section, we discuss how the model torus pa- 
rameters are constrained by our data set. As anticipated, we 
will consider all the 137 solutions at the l<x level. 

First of all, the torus is seen almost edge-on in all the so- 
lutions (i.e., the line of sight always intercepts the obscuring 
material), in agreement with the Type 2 X-ray and optical clas- 
sification of our sources. 

The torus model parameters which are left free to vary 
within the grid of models are: the torus opening angle 0, the 
slope a of the dust density profile and the optical depth T9.7. 

By converting the torus opening angle into a covering factor 
(CF) representing the fraction of solid angle covered by the 
dusty material, we find that solutions with high and low CF are 
possible, with a slight preference towards tori with large CF; 
the mean CF value is 0.65 (lcr~0.25), corresponding to a torus 
opening angle of ~1 10°. 



ties at 5100 A were derived from optical spectroscopy and the 
CF values from the ratios between the 6.7yum and the 5100 A 
luminosities for a sample of Type 1 quasars spanning five or- 
ders of magnitude in optical luminosity. In our approach, L5100 
is estimated, for each solution, from the input accretion-disk 
spectrum (see Sec. 14. Il l, once the normalization is found (see 
Sec. 16. U . The average value of CF an d L5100 for our sampl e 
lies, within l<x, on the relation found bv lMaiolino et al.1 J2007h . 
Unfortunately, given the limited range of L5100, we cannot in- 
vestigate the validity of the CF vs. L5100 relation over the range 
probed by Maiolino et al. (2007). 

Regarding the density profile, about 65% of the solutions 
have a — 0, while a = -0.5 and a = -1.0 represent 20% and 
15% of the solutions, respectively. This is reflected also in the 
16 best-fitting solutions, where only 2 sources (Mrk 0509#1 
and PKS 0537#111) are fitted with a = -0.5, one with a = 
— 1.0 (PKS 0537#43), and the remaining ones with a = (see 
TableE). 

Solutions with a moderate optical depth 79,7 are favoured 
by the SED-fitting analysis. As shown in ¥'\g.^\(bottom panel), 
there is a small number of solutions with high optical depths 
while the majority of the solutions (-80%) are characterized 
by 'moderate' T9.7 (r 9 .7<3) and 50% by low T9.7 (T9.7<1). The 
median value for T9.7 is 2. 

The finding of a preferred range of optical depths by 
the SED-fitting, even with relatively sparse photometric data, 
comes from the overall shape of the NIR/MIR continuum. In 
fact, once the stellar component is determined by the opti- 
cal/NIR data, the slope of the torus component is directly linked 
to the amount of absorption (i.e., to the optical depth) and is 
relatively well constrained by the available data. 

As shown in Fig.[3^,b, for the very low values of the opti- 
cal depth T9 7 (T9 j=0.1), the F06 model predicts spectra with a 
weak emission line at 9.7 //m. 

We clearly find that the optical depth T9.7 and the density 
profile a are not independent parameters, since low optical 
depth solutions mostly occur with flat density profile (a-Q). 
This is shown in Fig.\5\(top panel), where the fraction of solu- 
tions with a-Q is reported as a function of t 9 . 7 . For the assumed 
flared disk geometry, at high optical depth a flat density profile 
produce too much IR emission due to the large amount of dust 
at high radii. Thus the two best-fitting solutions with the highest 
T9 7 found (r 9 .7=6 and 10 for Mrk 509#01 and PKS 0537#111 
respectively, see Fig. [5] and Table [3) have a density profile de- 
creasing with the distance from the central BH (a=-0.5). 

Recalling that in our procedure different lines of sight are 
equivalent (having assumed y — in the radial density profile, 
Eq. 1), we can convert the optical depths to column densities 
N# (adopting a Galactic dust-to-gas ratio) to be compared with 
the Nh derived from the X-ray observations. Despite the un- 
certainties affecting the derivation of N# from the dust optical 
depths (i.e., dust and gas spatial distributions could be differ- 
ent), as well as tho se affecting the N# values from X-rays (see 



Pero la et al.1 12004 for details and Table 1), the two indepen- 



According to Maiol ino et al.l (120071) . the covering factor of 
the circum-nuclear dust decreases for increasing optical lumi- 
nosity at 5100 A (L5100) . This relation i s expl ained in terms 
of a "receding-torus". In iMaiolino et al.l (|2007) the luminosi- 



dent estimates give a consistent picture for the majority of the 
sources, once the lcr uncertainties, derived from the SED and 
X-ray fitting procedure, are taken into account. By excluding 
the two sources without a measured N# from the X-rays anal- 
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Fig. 5. Top panel: the fraction of solutions with flat density 
profile (a = 0) as a function of T9.7. Bottom panel: the T9.7 dis- 
tribution. The hatched and the empty distributions represent the 
best-fitting (16) and all the solutions (137) at l<x, respectively. 



ysis (see TableQ]), the median values for N# are ~7xl0 22 cm" 2 
and ~5.5xl0 22 cm" 2 , from the X-rays analysis and the dust 
optical depths, respectively. Therefore, the SED-fitting method 
confirms the X-ray classification of the sources as moderately 
obscured Compton-thin AGN. 

Two sources have a significantly different N# (by an order 
of magnitude) derived from the two methods, PKS 0537#1 1 1 
and Mrk 509#01. These objects are those characterized by the 
highest optical depths (T9.7 = 6, 10, which are converted into 
Nh~5.3x10 23 cm" 2 and N//~8.7xl0 23 cm" 2 , respectively). 



Since Mrk 509#01 has only an upper limit for the N# inferred 
from X-ray analysis (< 1.1 xlO 22 cm" 2 , see Table [TJ, the ob- 
served discrepancy for this object might be explained if the 
source is Compton-thick (Nh>10 24 cm" 2 ) and the observed 
X-ray spectrum is due to a reflection component. However, we 
cannot draw any firm conclusion on this issue. 

5.2. Host galaxy parameters 

In the spectral procedure, the host galaxy accounts for the 
optical/near-IR photometric data points, where the AGN contri- 
bution, given the obscured nature of our sources, is presumably 
low. We use the SSP spectra weighted by a Schmidt-like law 
of star formation (see Sec. l4.2l >. The extinction E(B-V) and r s f 
are free parameters. Once the best-fitting values for these two 
parameters are found, the stellar mass (obtained by integrating 
the star formation history over the galaxy age and subtracting 
the fraction due to mass loss during stellar evolution, ~ 30%, 
from it) is estimated from the SED normalization. At the end of 
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Fig. 6. 'Model' as a function of 'observed' bolometric lumi- 
nosities for the sample of luminous obscured quasars at the lcr 
significance level (137 solutions). Squares around the filled cir- 
cles represent the solutions with high optical depth (tc,^,,, > 3). 
The 'observed' bolometric luminosities are reported with no 
corrections applied. The solid line represents the identity re- 
lation. The dotted and dashed lines represent the predictions 
from the Pier & Rrolik (1992) model for 4 different configura- 
tions for Type 2 sources as a function of the optical depths t z 
and t r (where t z and r r are the optical depths at 9.7/im along 
the vertical- with respect the torus equatorial plane - and the 
radial direction, respectively) and viewing angle 9 (with re- 
spect to the equatorial plane). Red dotted lines: thinner model 
(t z = r r = 0.1, labeled with the letter t) with 6 = 0,50°. Blue 
dashed lines: thicker model (t z =1, T r =10, labeled with the let- 
ter T) with 6* = 0, 50°. At constant optical depth, the configura- 
tions with smaller viewing angle (closer to the equatorial torus 
plane) predict lower 'observed' luminosities. 



the SED fitting procedure, stellar masses are well constrained, 
for a given pair of E(B-V) and r S f , with a typical 1 <x uncertainty 
for the normalization of ~ 20%. 

We note that all but one of the stellar masses derived in this 
work are within 30% from the values estimated by Pozzi et al. 
(2007; eight sources in common), where the same data were 
used but different stellar libraries and a simpler approach was 
adopted (see Sect. 4.1 and 5.2 of Pozzi et al. 2007). 

The stellar masses found are in the range 4xl0 10 up to 
5xl0 12 M with three very massive galaxies (> 10 12 M Q , see 
Table [3), implying that our obscured AGN are hosted by mas- 
sive galaxies at hig h redshi f t. As s aid in Sec. 14. 21 the masses are 
obtained using a Isalpeterl d 19551) initial mass function (IMF) 
with mass in t he range ( 0.15-120 M B ). The assumption of a 
Chabrier (2003) IMF (see lRenzinill2006l) would produce a fac- 
tor of ~ 1.7 lower stellar masses. 

In Table [3] the best-fitting value for the free host galaxy 
parameters (r s f and E(B-V)) and the stellar masses are reported 
for each source. 
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6. Black hole physical properties 

6.1. Black hole accretion luminosities 

The accretion-disk luminosity L acc is a direct output of the fit- 
ting procedure and is obtained by integrating the code input 
energy spectrum once the best-fitting torus components and its 
normalization are found (see Sect. 4.3). 

The input spectrum is defined in the 1CT 3 to 20 fim regime. 
Although this wavelength range provides the largest contribu- 
tion to the nuclear AGN luminosity, we have also included in 
the L acc computation the hard X-ray luminosity (L i.25-500fe v)- 
This luminosity is estimated from the de-absorbed, k-corrected 
La-Wkev luminosity, assumping a photon index T=1.9 (typical 
of AGN emiss ion) and an exponential cut-off at 200 keV (e.g., 
Gilli et alJ2007h . A different choice for the energy cut-off (e.g., 
at 100 keV) would produce a difference by « 20% in the to- 



tal X -ray luminosity for F= 1 .9 sources (see IVasudevan et al 
20091) . 



We note that dust grains are almost transparent to hard 
X-ray photons, therefore the output of the code is not affected 
by the fact that the accretion-disk model spectrum does not ex- 
tend above soft-X-ray energies. 

In Table [3] L hoh along with L acc and Li.25-500teV. are re- 
ported. Lboi extends over two orders of magnitudes (10 44 - 
10 46 erg s _1 ), with the hard-X-ray luminosities (1.25-500 keV) 
contributing in the range 5-50% of the AGN power (see 
Table [3]). The two sources with the highest optical depths 
(i'9.7 jUm =6,10) are among the sources with the smallest hard- 
X-ray fraction (Mrk 509#01 and PKS 0537#111). In Table [3] 
we report also the range of bolometric luminosities as obtained 
by considering the full set of lcr solutions. The uncertainties 
are, on average, of the order of 0.2 dex, but vary significantly 
from source to source, ranging from about 5% to about a factor 
3 (see also Fig. [7]). 

We compare the computed bolometric luminosities with the 
luminosities derived by integrating the torus best-fitting tem- 
plates from 0.1-1000/vm (plus adding the hard X-ray luminos- 
ity for self-consistency). The two methods assume the same 
torus SED, hence the comparison can give important informa- 
tion on the systematics affecting the estimates of Lt, i derived 
by integrating the observed SED, which is the method widely 
used in literature. We refer to the first measures as the 'model' 
luminosities and to the latter as the 'observed' luminosities. 

The 'observed' Lboi (see Fig.|6]l are lower (up to an order of 
magnitude) than the 'model' ones for all the solutions; the me- 
dian value of the ratio is R~2 (the solid line in Fig.|6]represents 
the identity relation). An under-estimate of the luminosity in 
Type 2 sources is expected by to rus models (e.g. JPier & Krolikl 
1993c iGranato & Danesdll994l) : here, we quantify this effect 
and provide an empirical factor to correct the 'observed' lumi- 
nosities, at least for this class of sources. 

We underline how the observed discrepancy does not de- 
pend on the lack of observations at far-IR wavelengths. In 
fact, the two methods assume the same torus SED for self- 
consistency (i.e., the integrated torus SED to estimate the ob- 
served IR luminosity is the output of the code); under this hy- 
pothesis, an over(under)-estimate on one luminosity would in- 



troduce the same effect on the other. As a result, a poor sam- 
pling in the far-IR would have the same impact on both (i.e. 
'observed' and 'model' ) luminosities. Our analysis takes into 
account this uncertainty by considering all the solutions (i.e., 
all torus models) at the lcr confidence level. By means of this 
procedure, a broad range of model SED is associated to each 
source (on average, eight solutions; see Sec. [5]), characterized 
by different emission in the mid/far-IR region, as a result of 
different torus geometry and absorption properties (see Fig. [3). 

As explained in Pier & Krolik (1992), the low values of 
the 'observed' Lboi depend on a combination of three factors: 
the torus opening angle (geometrical factor), the observer 
viewing angle 6 and the torus optical thickness Tgj^. By erro- 
neously assuming isotropic torus emission (as done to compute 
the 'observed' L/, /), the primary flux which does not intercept 
the obscuring material would not be included in the luminos- 
ity budget; moreover, as the thickness of the torus increases, 
more and more primary high-energy photons entering the torus 
are absorbed by the dust grains and re-emitted isotropically 
(hence also in directions escaping the torus it self). This ef- 
fect is explained by the dust self-absorption, i.e. thermal dust 
emission absorbed by the dust itself. For high optical depth, the 
outer edges of the torus absorb the IR photons coming from the 
warmer dust at smaller radii and re-emit them isotropically, i.e. 
also in directions outside the line of sight. To better visualize 
this effect, we report in Fig. [6] the 'observed' versus 'model' 
luminosities, as predicted by Pier & Krolik (1992), as a func- 
tion of the viewing angle 6 and the torus optical thickness Tg.j Mm 
for 4 sets of Type 2 configurations (as described in Fig. [(^cap- 
tion). Although there are some slight differences between the 
F06 model (adopted here) and the Pier & Krolik (1992) torus 
model (where the optical depth varies independently along the 
radial and the vertical axis), optically thinner models show less 
anisotropy (red dotted lines in Fig. [6] corresponding to two dif- 
ferent viewing angles), than higher Tg.T^ m models (blue dashed 
lines in Fig. [6] corresponding to the same viewing angles con- 
sidered for the thinner model). The cold outer edges of the 
thicker models, in fact, radiate little and block the light com- 
ing from the inner torus radii. 

To investigate these issues further, we apply a 'conserva- 
tive' correction to our 'observed' luminosities, meant to correct 
only for the geometrical factor; in other words, we divided each 
'observed' luminosity by the corresponding covering factor CF 
(-0.58 for -100° and -0.88 for -140°). Although this correc- 
tion increases the 'observed' luminosities, the 'model' ones are 
still higher (R~l.6); the remaining discrepancy is mostly found 
for solutions with high optical depth, as expected (R ~ 5 for 
models with T9 7>3; see Fig.|6l where the squares mark the 52 
solutions with t 9 7>3). 

An independe nt and consistent analysis was done also by 



Pozz i et alJ (120071 see their § 5.1) where a first-order correc- 
tion of ~2 to the 'observed' luminosities was estimated, ac- 
counting for geometrical and anisotropy effects; in that work, 
however, the correction was estimated usin g the ratio of ob- 
scured/unobscured quasars according to the iGilli et a l. (2007) 
AGN synthesis models of the X-ray background and the differ- 
ent shape of Type 2 vs. Type 1 quasar SEDs as a function of the 
column density. In lPozzi et al J ( 12007b . the SED fitting was done 
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Fig. 7. 2-10 keV bolometric corrections as a function of the 
'model' bolometric luminosities (filled circles). Filled circles 
inside empty red squares represent the sources with a spectro- 
scopic redshift. The r ed solid and do t -dashe d lines represent the 
predictions from the iMarconi et al. I d2004 relation and its lcr 
dispersions. Also the expectations from iHopkins et al.l (12007b 
are reporte d as empty blue trian gles. The red dotted line rep- 
resents the IMarconi et al.l (120041) expectations at 5 <x from the 
best-fitting relation. 



using the lSilva et al.l ( 120041) AGN templates. Since the template 
choice was based on the X-ray N# (and not on the N# resulting 
from the torus modelling as in the present analysis), the correc- 
tion corresponding to the thicker models (Nh 1 10 24 crrT 2 ) were 
not included since no Compton-thick objects were revealed in 
X-rays. 

6.2. Hard X-ray bolometric corrections 

In Fig.Qthe bolometric-to-X-ray luminosity ratio (fe-iofev) as 
a function of Lj ; is shown: for the bolometric luminosities we 
assume the model ones. The error bars on &2-iofev are derived 
from the lcr dispersion on L\,„i. 

A large spread in the fo-iofev is found ( 6;fc2-iofa?y ~ 80), as 



pointed out also by the pioneering work of Elvis et al .1 (11994) 
on Type 1 QSOs, due to the large dispersion in the AGN spec- 
tral shape. Our median value (&2-iofev ~ 20, estimated on the 
137 solutions), is marginally consistent with the mean value of 
Elvis et al.l d 1994b . ~25 after removing the IR contribution (in 
order to not double-count the fraction of the nuclear emission 
absorbed by the circumnuclear dusty material seen almost face- 
on). We confirm the trend of higher A^-io^y for objects with 



higher bolometric luminosities as predicted by Marc oni et al 



(2004) (red solid line in Fig. [2 red dot-dashed lines represent- 
ing the lcr model dispersion), but our kz-\okeV values are sig- 
nificantly lower (at least a factor 2 in normalization). They de- 
rive k2-\ok e v by constructing an AGN reference template taking 
into account h ow the spectral inde x a ox varies as a function of 
the luminosity dVignali et al]2003l) . Predictions consistent with 
Marconi et alj ( 120041) were obtained more recently by Hopkins 



et al. (2007, blue triangles in Fig. [7J, cons idering the most re- 
cent determ ination of SEP tem plates (i.e. lRichards et al . 2006) 
and a ox (le. lSteffen etaDl2006l) . 

Our low values for k^-iokeV a re consistent wit h our previous 
analysis (median fo-iofev ~ 25. iPozzi et al.ll2007l) based on a 
differ ent method and on different AGN templates ( Silva et aD 
2004) and with other estimates found in literature for hard 



X-ray selected samples. Kuraszkie wicz et al .1 (120031) . consid 



,43 



ering a sample of X-ray selected luminous AGN (10 
< 10 46 erg sec ) found a median &2-iofev 

of 18; 



Ball o et alJ (120071) . analysing a sample of low-luminosity AGN 
(10 42 < La-iokev < 10 43 6 erg sec -1 ), found a median ^-lOfev 
of 12. Low bolometric-to-X-ray ratios, consistent with our es- 
timate (med i an fe- iofey ~ 25, lcr=53) were found recently by 
Lusso et al. (2009), where the statistical properties of a large 
(and complete) sample of 545 X-ra y selected Type 1 Q SO 
from the XMM-COSMOS survey faasinger et all 120071) are 
presented. 

The lower bolometric-to-X-ray luminosity ratios found in 
the present work ( and in the abo v e men tione d samples), in 
compa rison to the Marconi et al.l d2004l) and Hopkins et al 



(2007) predictions, are probably caused by a selection bias, 
since our sample (and most of the above cited ones) are hard- 
X-ray selected samples (i.e., sources with high X-ray luminos- 
ity are favored). Moreover, as discussed in §|2j our sources are 
among the most extreme X-ray sources, being characterized by 
red optical-to-NIR colours (R-K s ; 5) and high X-ray-to-optical 
ratio (X/O; 1). Our selection is likely the origin of the large de- 
viation (at about the 5<x level) f or a large fraction of th e present 
sample (see Fig. |7J from the IMarconi et al.l (|2004) relation. 
To further explore this issue, a larger (and complete) sample 
of X-ray sources (with optical identification up to the faintest 
X-ray fluxes) is needed, in order to correct for the selection bias 
and to derive the properties of the parent AGN population (see 
Lusso et alj|2009b . 



6.3. Black hole masses 

The BH masses are not a direct output of the best-fitting pro- 
cedure and cannot be derived using 'standard methods' (i.e. 
galaxy stellar kinematics, nuclear gas motions, reverberation). 
We estimate th em indirectly using the M buhe-MsH relation de- 
rived locally bv lMarconi & Hundd2003l) . by assuming as M/,„/ ge 
the stellar mass derived from our best-fitting procedure. The 
main uncertainties affecting these estimates derive from the ex- 
trapolation of the local relation to higher z, where the behaviour 
of this relation i s still a matter of debate (see discussion in 
Pozz i et al. 1 120071) . As far as the stellar masses are concerned, 
they are quite well constrained by the SED-fitting procedure 
inside the pre-constructed grid of galaxy models (see Sec. 15. 2b . 

The inferred black hole masses are typically in the range 
10 8 - 1O 9 M , with three sources (PKS 0537#43, GD 158#19, 
Abell 2690#29) with higher masses (M BH ~10 9 - 5 - 1O 1O O M ). 

The r ange of BH masses i s con sistent with the values re- 
ported by iMcLure & Dunlopl d2004 for the SPSS quasars in 
the same redshift interval (0.9 t z 1 2. 1 , see also Sh en et al . 2008, 
where new BH masses are derived). 
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l2010h . the Eddington ratios Xem would decrease further by the 
same factor. 

Finally, in Fig. [9] the bolometric corrections fe-iofev 
are plotted against the E ddington ratios Aem (following 
Vasudevan & Fa bian] l2009l) . The sources of the present work 



are reported as black filled circle (the error bars representing 
the lcr confidence interval derived from t he uncertainties on 
Lboi)- Along with our data we show the IVasudevan & Fabian] 
(2009) results, where simultaneous optical, UV and X-ray ob 



servati ons are included for the majority of the iPeterson et al 



0.0 0.5 1.0 1.5 
z 



2.0 



(2004) reverberation mapped sample of AGN (blue empty cir- 
cles). Our findings are in fairly good agr eement with the trend 
of inc reasing ki-mkev for increasing Aem- 1 Vasudevan & Fa bian 
(2009) interpret the observed trend as due to different black 
hole SED shape as a function of the Eddington ratio, with the 
high and low Eddington ratios corresponding to different frac- 
tions of the ionizing UV luminosity. A similar trend was re- 
cently found for a sample of 63 Type 1 and Type 2 AGN de 



Fig. 8. Aem as a function of z. Black circles: sources of the 
present sample. Red squares as in Fig. 7. The error bars rep- 
resent the lcr uncertainties on Lj ; (as derived from the x 1 
analysis). Small grey cro sses: sample of SDSS quasars from 
McLure & Dunlopl (2004), with the median values (and asso- 
ciated uncertainties assuming the normalized median absolute 
values) of Aem within Az = 0.1 bins shown as red filled circles. 
Blue triangles: median values (and associated uncertainties as- 
suming the normalized median absolute values) of Azdd within 
Az = .2 bins from the COSMOS Type 1 sample o flLusso et al 
d2009h . 



tected in the Swift/BAT 9-months catalog (see Vasudevan et al 
2009). At variance with the assumption in lVasudevan & Fabian 



(2009), where the bolometric luminosities were derived by 
integration over the observed optical/UV/X-ray SED, in this 
work the authors consider the reprocessed IR emission, repro- 
duced by the empirical SEDs of Silva et al. (2004), as a proxy 
of the intrinsic AGN bolometric luminosity, as firstly suggested 
by Pozzi et al. (2007). 

The dependence of bolometric corrections on Eddington 
ratios is expected by accretion-disk models, which predict an 
increasing hard X-ray bolometric corrections at increasing ac- 
cretion rates (e.g. Witt et all 1997 ). Recently also Bianchi et al 



In Fig. [8] the Eddington ratios Aem< defined as Aem - 
L b oilL E dd (with L Ed d = 1.38xlO 38 M B// /M ), are reported 
as a function of redshift. The values are compared with 
those of the whole SPSS quasar sample (small g rey crosses 
McLure & Dunlopl 2004) and those obtained by Lusso et al 



(2009), studying a large (156 sources) sample of Type 1 X-ray 
AGN from the XMM-Newton archive suggest that the bolo- 
metric correction must depend on Eddington ratio in order to 
allow the intrinsic power of AGN to scale linearly with black 
hole masses. 



2009) for the sub-sample of 150 X-ray selected Type 1 AGN 7 - Summary 



in COSMOS with an accurate black hole mass determination 
(blue triangles). 

The Aem values of the present work cover slightly more 
than an order of magnitude (0.01-0.3), with a median value of 
AEdd~0-08 (estimated considering all the 137 solutions at lcr 
level, see Sec. f5J. The derived values are within the 3cr confi- 
dence interval of the SDSS quasar Aem distribution, charac- 
terized by a median value of ~0.3 and with a dispersion of 
0.35 dex at the same redshift interval sampled by our sources. 
However, almost all our data points lie towards the low Aem 
tail of the SDSS distribution (see Fig. [8), suggesting that X-ray 
selection is biased towards slightly lower Aem man optical se- 
lection. Our data are fully consistent with the results obtained 
from a much larger sample of X-ray selected Type 1 AGN in 
the COSMOS field dLusso et al.ll2009h . 

The results are robust against the uncertainties on the ex- 
trapolation, discussed above, of the local Mbuige - Mbh rela- 
tion, at the redshift of our sample. In fact, allowing for positive 
evoluti on with redshift of the Mbh I Mbuige ratio by a factor of 
2 (e.g jHopkins et al.ll2006T: IShields et atl l2006: Mer loni et al 



We analyzed the SEDs of a sample of 16 obscured quasars se- 
lected in the hard X-ray band. Spitzer mid/far-IR photometry 
(IRAC and MIPS), along with the data available in the liter- 
ature, is modeled using a multi-component model, where the 
AGN re-processed emission is reproduced in the context of a 
flared disk model, as described by F06. Within the context of 
a flared disk torus model, the uncertainty in and degeneracy 
between the various derived parameters are accounted for by 
including all solutions within lcr of the best- fit in the subse- 
quent analysis. 

The main results are summarized below: 

• All the 16 quasars are detected up to 8 /mi and all, but 
two sources, are detected at 24/mi with flux densities in the 
range 100-7000 /vJy at the 5<x level. The two most luminous 
sources of the sample are detected also at 70 and 160 fim. 

• The observed broad-band spectral energy distributions are 
well reproduced by a multi-component model comprising 
a stellar, an AGN and a starburst components (when far-IR 
detections are available). The AGN component, modelled 



PKS 0312#36 
PKS 0537#91 
PKS 0537#54 
PKS 0537#111 
Abell 2690#29 
PKS 0312#45 
BPM 16274*69 



9.6/5 
2.9/5 
2.9/5 
11.4/5 
0.5/2 
0.5/5 
0.5/5 



0.0 
0.0 
0.0 
-0.5 
0.0 
0.0 
0.0 



0.1 
3.0 
0.1 
10.0 
1.0 
0.1 
0.6 



140 
140 
140 
140 
140 
100 
140 



0.85 
0.25 
0.05 
0.05 
0.15 
0.25 
0.05 



0.79 
0.3 
0.17 
0.31 
0.55 
0.48 
0.25 



3.0 
66.0 

5.6 
103.1 
168.8 
24.25 

11.4 



2.5 
23.7 

4.6 
13.4 
94.7 
11.3 

7.3 



0.7 
8.1 
2.0 
1.7 
8.4 
6.2 
2.4 



2.7 
31.4 
7.8 
6.6 
32.6 
24.0 
9.3 



5.7 [5.7-8.1] 
98.3 [40.2-326] 
13.3 [13.3-25.8] 
110 [34.1-110] 
201 [135-264] 
48.3 [37.4-239] 
20.7 [18.2-56.6] 



8.2 
12.1 
6.7 
64.4 
23.9 
7.8 
8.6 



0.43 
1.48 
1.72 
4.61 
17.9 
3.61 
4.68 



0.11 
0.35 
0.40 
1.03 
3.80 
0.81 
1.05 



0.043 
0.224 
0.026 
0.084 
0.042 
0.047 
0.016 



(1) source name 

(2) best-fitting minimum^- 2 0f 2 ,„) and number of parameters to be fitted (degrees of freedom); 

(3) , (4), (5) best-fitting torus parameters (a: exponent of the power law parameterizing the density profile; r 91llm : 9.7 /zm optical depth, 0: torus opening angle). The ratio 
parameter R ma . x /R„,i„ is frozen to 30; the density parameter y is frozen to 0; 

(6), (7) best-fitting stellar parameters ij s f. duration of the exponential decay of the burst in units of the oldest SSP; E(B - V): extinction); 

(8) accretion-disk model luminosity, L acc (from soft X-ray to IR frequencies) which represents the torus model input luminosity, in units of 10 44 erg s" 1 ; 

(9) integrated (1-1000/im) torus luminosity in units of 10 44 erg s -1 (not corrected, see l6.lt ; 

(10) absorption-corrected 2-10 keV luminosity in units of 10 44 erg s _1 ; 

(11) hard-X-ray (1.25-500 keV) luminosity in units of 10 44 erg s -1 (derived from the 2-10 keV luminosity, see !6.2b : 

(12) bolometric AGN luminosity (L acc + Li.25-5oofey) m units of 10 44 erg s -1 ; the lcr range derived from the SED-fitting analysis is reported; 

(13) 2-10 keV bolometric correction (L bo il L 2 -mev)\ 

(14) galaxy mass in units of 10 u M o ; 

(15) black hole masses in units of 10 9 M o (estimated from lMarconi & Huntll2003l relation); 

(16) Eddington ratios (L bo ,/L Edd ). 
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Fig. 9. k2-\okeV as a function of Asdd- Black filled symbols: 
sources of the present sample. Red squares as in Fig. 7. The 
error bars represent the lcr uncertainties on Lt, i (as derived 
from the^ 2 analysis) which affect both th e kx-wkev and the Azdd 
value s. Blue open symbols: sources from lVasudevan & Fa bian 
(2009). 



with the F06 radiative transfer code, accounts for the X-ray 
emission and for a fraction of the IR emission, mainly due 
to reprocessed emission from the putative dusty torus sur- 
rounding the central black hole. 

Solutions with a moderate optical depth T9.7 are favoured 
by the SED-fitting, with the majority of the sources hav- 
ing moderate optical depths (j^j^-^i). The derived gas 
column densities (N#) are consistent, for most of the 
sources, with the values estimated from the X-ray analy- 
sis, both indicating that the sources are Compton-thin AGN 
(N H ~10 22 - 3xl0 23 cirr 2 ). 

The model nuclear bolometric luminosities are in the range 
5xl0 44 - 4xl0 46 erg s _1 . By comparing these values with 
those obtained by the integration of the nuclear observed 
SED, we conclude that the latter under-estimate the bolo- 
metric luminosities by a factor of 2. The difference may be 
explained by anisotropic torus emission and the effect of 
the torus optical depth (e.g., Pier & Krolik 1992). 
From the model nuclear SEDs, we estimate the bolometric- 
to-X-ray corrections (fo-iofev)- The median is ~20 

(6;A:2-io*:cy~80). The value is smaller than assumed by 
some models of BH evolution (fe-iofevMO at the median 
luminosity of our sample). The discrepancy is significant at 
5<x level at low bolometric luminosity. 
By assuming the local Mt, u i ge - Mbh relation, we estimate 
A-Edd with a median value of 0.08 (0.01 tA,Edd~ 0.3). The 
whole SDSS quasar sample, at the same redshift interval 
sampled by our objects, is characterized by a median value 
of 0.3. Our data are within the 3cr confidence level of the 
optically selected quasar distribution. However, almost all 
our sources lie towards the low Xem tail of the SDSS dis- 



tribution, suggesting that our X-ray selection is biased to- 
wards lower Eddington efficiencies than optical selection. 
• The data are consistent with the correlation recently sug- 
gested by Vasudevan & Fabian (2007, 2009) between 
ki-wkev and A E dd, where low bolometric corrections are 
found at low Eddington ratios. 
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